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Abstract 



We study the matter effects for solar neutrino oscillations in a general scheme, without any 
p£| ■ constraint on the number of sterile neutrinos and the mixing matrix elements, only assuming a re- 

alistic hierarchy of neutrino squared-mass differences in which the smallest squared-mass difference 



Oh 1 



Qh| is effective in solar neutrino oscillations. The validity of the analytic results are illustrated with a 



numerical solution of the evolution equation in the simplest case of four-neutrino mixing with the 
realistic matter density profile inside the Sun. 



j> . Neutrino physics [1-3] is one of the most active fields in particle physics. The standard scenario [3-5] 
! in neutrino oscillation phenomenology is three neutrino mixing with a squared-mass hierarchy: 

Am^ OL ~ 8xl(T 5 eV 2 , Am\ TM ~ 2.5xl(T 3 eV 2 . (1) 

The existence of a much larger squared-mass difference ( > 0.1 eV 2 ), as that indicated by the LSND [6] 
Q\ ', — > v e oscillation signal, would require the existence of sterile neutrinos in addition to the three active 
^ \ flavors. However, the LSND signal is currently disfavored by the negative results of the KARMEN [7] 
^ ■ and MiniBooNE [10] experiments. Another indication comes from the anomalous ratio of measured 
^ ■ and predicted 71 Ge observed in the Gallium radioactive source experiments GALLEX [8] and SAGE [9] 
and the MiniBooNE [10] low-energy anomaly, which can be explained by Short Baseline (SBL) electron 
neutrino disappearance [11]. These possible active-sterile transitions can be tested by studying their 
effects in solar neutrino oscillations. 

The matter effects in solar neutrino active-sterile oscillations were studied in Ref. [12] in a four neutrino 
scheme, but the SBL effects were neglected. Thus a combined analysis of the two anomalies above with 
the solar neutrino data needs the derivation of the matter effects in a more generic scheme [13]. 
To start, we consider a generic scheme of mixing of three active neutrino fields (v e , u^, u T ) and N s sterile 
neutrino fields mixed as v aL = J2k=i Uak^kL (oc = e, /i, r, si, . . . , sn s ) , with N = 3 + N s . Solar neutrinos 
are described by the states = Y^ a i'ai^Wa) , where x is the distance from the production point 

with Va(0) = 5 ae , and the amplitudes are normalized as J2 a li'ai^l 2 — 1 • The evolution of the flavor 
transition amplitudes ip a (x) is given by the MSW equation (see Ref. 



i±* = ±(UM>tf + A)1, (2) 
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where E is the neutrino energy and 



(ip e ,ij fM ij T ,ij sl ,...,ip SNs )' 



M 2 = diag(0, Am 2 1? Am 2 31 , A 



m 



41' 



Am 2 m ) 



A = diag(A C c + ^nc, ^nc, -4nc, 0, . . .) 



(3) 
(4) 
(5) 



with Amjy = m 2 — m 2 the squared-mass differences and A cc = 2EVqc , A^c — 2EVnc ■ ^cc an d Vnc 
are the charged- current and neutral-current potentials. 
For solar densities we have 



.4 



cc 



1-4 



NC 



Am^ < | Ami, | for k > 3 . 



'-fci 



(6) 



It is useful to work in the vacuum mass basis \l/ v = Ufi[ ■> ■ ■ ■ > ? /'^) T = > which satisfies the evolution 
equation: 

,i_v|/V = (m 2 + C/UC/) M> v . (7) 



dx 



2E 



The inequality in Eq.flS]) imply that the evolution of the amplitudes ipj, . . . ,ipjf is decoupled from the 



others. Then we have 



^)-^ V (0)exp(-z^|^ 



for k > 3 



with ip^(0) = U* k . The first two amplitudes are coupled by the matter effects: 

d fipY\ 1 /- Am\ x + A cos 2£ A sin 2£ \ (ipj 



dx 



AE 



A sin 2£ 



Ami — A cos 2£/ 



(9) 



with the definitions of tan 2^ = Y/X , A = A cc \/^ 2 + Y 2 , and 



X = \U el \ 2 - \U e2 \ 2 + Rvc (l^iT- l^l 2 ) , Y = 2 



U* 1 U e 2 + RnC 2^ Ual&a2 



(10) 



From the similarity of ([9]) and the corresponding equation for v^-v^ or v e -v T two-neutrino mixing (see 
Ref. [2]), we obtain the averaged oscillation probability: 



P 



~ + Q-P12) cos 2(^00^2/' 



N 



COS 2 Xp COS 2 + \ U Pk\ 2 \Uek\ 2 , 



11^ 



with the mixing angles dp and Xp defined by \U f 



0i\ 



k=3 

cos 2 tip cos 2 X/3 , 1^/32 1 2 = sin^i?^ cos 2 X/3 , and 



sin 2 = X]fc=3 l^/3fc| 2 • The effective mixing angle in the production region can be written as "d° e = ~d e 
with tan2u; = (A sin 2£)/(Am 2 1 — y4°cos2£) . The crossing probability P\i is given by 

p _ ex P(-f^)-ex P (-|7R-fe) 

"12 — 7 n v{A - A R J , 

1-exp (-|7r-^J 

where 7r is the adiabaticity parameter at the resonance (A-r = Am^ cos 2£) : 

Am 2 sin 2 2£ 

7r 



-oj 



'121 



(13) 



2£cos2£|dlniV e /d£| R ' 

To check the validity of the results for the realistic BP04 Standard Solar Model density [14] , the results 
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Figure 1: Averaged probability of u e survival and v e — > v s transitions as functions of the neutrino 
energy E for the mixing matrix in (Till) calculated for the BP04 Standard Solar Model density [14]. The 
lines are obtained with the analytic expression in Eq. (1111) and the overlapping points are obtained with 
a numerical solution of the evolution equation. 



of the analytic expression (TTTj) and of the numerical solution of the evolution equation ([9]) are presented 
in Fig. CD The relevant mixing parameters are chosen as tan 2 d e = tan 2 i^sol — 0.4 and 



\U, 



<-■'-> I 



0.05, \U. 



e -l | 



0.05 



Wsl\ 



0.03 



\u, 



,s2 



0.06. 



(14) 



We considered only neutrino energies smaller than about 50 MeV, for which the inequality ([6]) and 
the approximation ([8]) are valid. In this range of energies, from Fig. [1] one can see that the analytic 
approximation is very accurate for solar neutrinos and can be safely used in the analysis of solar data [15]. 
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